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EXPERIMENTAL TECHNIQUES FOR EVALUATING STEADY-STATE JET 
ENGINE PERFORMANCE IN AN ALTITUDE FACILITY 
by John M. Smith, Chi Y. Young, and Robert J. Anti 
Lewis Research Center 

SUMMARY 

Airflow and thrust calibration tests were conducted on a J85 afterburning turbojet 
engine in an altitude facility of the NASA Lewis Research Center using a contoured bell- 
mouth inlet duct, four fixed -area water-cooled exhaust nozzles, and an accurately cali- 
brated thrust system. 

Accuracy in these calibrations was of prime concern in this program. The probable 
experimental error in airflow measurement was 0. 6 percent near the bellmouth throat 
design Mach number of 0.6. This probable error was based on individual and cumula- 
tive errors of instrumentation. The probable error in gross thrust measurement was 
approximately 0. 6 percent at the load cell full-scale reading. This was based on in- 
strumentation, data systems, and thrust load cell errors. 

A simplified method for calculating airflow rate is presented and validated by a 
more rigorous method when the inlet airflow conditions are favorable. 

The isentropic thrust coefficients and the gross thrust coefficients for the exhaust 
nozzles were experimentally determined with and without an engine in place. A J85 
turbojet engine was used for the hot runs . 

The values of the isentropic thrust coefficients from hot -gas data were on the aver- 
age within ±0. 59 percent of the cold-flow values using the cold-flow data as a base line. 


INTRODUCTION 

Flight flying test bed and wind-tunnel evaluation of supersonic exhaust nozzles oper- 
ating in the transonic flight regime are currently underway at Lewis. In order to deter- 
mine nozzle characteristics, the nozzle inlet total pressure, temperature, and gas flow 
rate must be known. Obtaining these measurements is sometimes difficult as in the case 



of the flying test bed. Therefore, these variables must be determined accurately in 
ground -test facilities. 

Engine thrust is one of the prime design variables and accurate calibrations and 
measurements are needed to determine the gross thrust and nozzle thrust coefficients. 

A prime requirement is the design of a thrust measuring system to evaluate gross thrust 
within the required accuracy. 

Airflow is also a prime engine variable. An airflow measuring station is invariably 
incorporated in any engine test setup. The flow calibration of the airflow measuring sta- 
tion can be determined using temperature and pressure measurements. 

Flow conditions at the inlet of the exhaust nozzle must be known in order to deter- 
mine the nozzle coefficients. The use of fixed-area water-cooled nozzles provides a 
means of accurately determining nozzle throat area by minimizing thermal expansion. 
The technique of using water-cooled nozzles is reported on in reference 1. Exhaust noz- 
zle thrust coefficients for both hot -gas and cold-flow rims have been calculated and com- 
pared using the cold flow data as a base line. This comparison is needed to show that 
the thrust measuring system designed for these tests could evaluate gross thrust within 
the desired accuracy for a turbojet engine, and that the thrust system could also be used 
for other engine and nozzle evaluation tests. Several nozzle evaluation test results are 
reported in references 2 to 4. 

Hot gas was provided by a J85 turbojet engine with an afterburner with exhaust gas 
temperatures up to 1944 K (3500° R). 

All tests were conducted in the NASA Lewis Research Center Propulsion Systems 
Laboratory . 


APPARATUS 

Figure 1 shows a schematic diagram of the contoured inlet air bellmouth that was 
used for the tests . The inlet of the bellmouth and a portion of the length was based on an 
ASME design. The remainder of the length was designed for sinusoidal acceleration. 

The inlet airflow measuring station (station 1) was designed for a throat Mach number of 

0 . 6 . 

Figure 2(a) shows a schematic drawing of the test installation of the duct used for 
the cold-flow tests. The dimensions of the straight cold-flow duct were approximately 
the same as a J85 turbojet engine. It was 40.64 centimeters (16 in.) at the inside diam- 
eter and 238.76 centimeters (94 in.) long. A screen was installed in the duct for use 
with any nozzle that had a throat diameter larger than that of the bellmouth straight sec- 
tion. This was to insure that the nozzle would choke instead of the bellmouth straight 
section. 
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Figure 1. - Inlet air bellmouth. 


After the cold -flow tests, a J85 turbojet engine was installed as a gas generator for 
the hot -gas tests. Figure 2(b) shows a schematic drawing of the test installation of a J85 
engine for the hot -gas tests. The engine was approximately 50.8 centimeters (20 in.) at 
the outside diameter and 276.86 centimeters (109 in.) in length including the afterburner. 
The engine has an eight-stage axial-flow compressor coupled directly to a two-stage tur- 
bine. The engine also has an annular primary combustion system, variable inlet guide 
vanes, controlled compressor interstage bleed, afterburner, and variable -area exhaust 
nozzle. The variable-area nozzle was replaced by the fixed-area water-cooled nozzles 
used for these calibration tests and those reported in reference 1. 

Figure 3(a) is a photograph showing a water-cooled nozzle with cooling lines at- 
tached, fastened to the afterburner of the J85 engine. A schematic of the interior con- 
tours of these nozzles is shown in figure 3(b). The nozzles were convergent-divergent 
types with a contoured convergent section and a throat radius of curvature equal to the 
throat diameter. The divergent portion was a 15° half-angle cone with an expansion area 
ratio of approximately 1.3. The throat areas were 694, 843, 1073, and 1186 square 
centimeters (107. 5, 130.7, 166.3, and 183. 8 in. ^), respectively. 
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(a) Test installation of duct for cold-flow tests. 



Figure 2. - Test installation for cold-flow and hot-gas tests. 
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FACILITY 


The test pregram was conducted in the Propulsion Systems Laboratory Altitude 
Chamber at the NASA Lewis Research Center. In this facility, altitude from sea level 
to 24 000 meters (80 000 ft) can be simulated with airflow rates from 218 to 23 kilograms 
per second (480 to 50 lb/sec), respectively. Pressure ratio can also be varied over a 
wide range. In this test program the inlet to exhaust pressure ratio covered a range 
from 1. 2:1 to 32:1. The inlet air bellmouth passed through a labyrinth seal in the cham- 
ber forward bulkhead and was directly connected to either the straight duct or the J85 
engine. This system was mounted on a thrust measuring stand as shown in figure 2. 


INSTRUMENTATION 

Instrumentation was provided to measure the necessary parameters to calibrate the 
airflow measuring station, evaluate the flow and thrust coefficients of the nozzles, and 
evaluate the experimental and theoretical terms in the thrust equation. Measurements 
were located in the inlet plenum P the airflow measuring station (1), the nozzle inlet 
(7), and the altitude chamber (0) as shown in figure 2. Details of the locations of the 
pressure and temperature instrumentation at each station and on the bellmouth surface 
are shown in figure 4 as viewed from upstream. Other pressure probes located at the 
back of the bellmouth slightly ahead of the labyrinth seal and in the altitude chamber are 
not shown. 

The data were recorded on the Lewis central automatic digital data encoder 
(CADDE) . Temperatures were recorded as voltage outputs on the automatic voltage digi- 
tizer (AVD) system, and pressures were recorded on the digital automatic multiple pres- 
sure recording system (DAMPR) . The data were then processed in a digital computer. 

A detailed description of the CADDE system is given in reference 5. 

The total pressure probes used on the pressure rakes were single tubes with a 40° 
internal chamber. Information on the design of these probes can be found in reference 6. 
Total temperatures at stations 1 and 7 were corrected using recovery factors that were 
function of Mach number. These recovery factors were obtained from data in refer- 
ence 7. 
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(a) Inlet plenum (sta. PI cross-rake, 30.48 cm (12 in.) (b) Bellmouth static pressure, 

upstream of bellmouth. 


Total pressure and 
temperature rake^ 



Water-cooled 
total pressure rake^ 



(c) Station 1 - downstream end of bellmouth throat (01 Nozzle inlet (station 7). 

straight section. 

CD-I 1069 -28 


Figure 4. - Bellmouth and duct instrumentation stations as viewed from upstream. 


THRUST SYSTEM 

One of the primary purposes of these tests was to determine the accuracy of the 
thrust measuring system in evaluating gross thrust for the turbojet engine and whether 
this thrust system could be used in other engine and nozzle tests . 
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As can be seen in figure 2, the test bed supporting either the straight duct or the J85 
engine was suspended by flexure plates from facility supports. 

A Baldwin Lima Hamilton double bridge load cell with a 22 240-N (5000-lbf) range 
was used to measure thrust. It was water cooled to minimize drift from heating in the 
test cell. The load cell was calibrated with a Morehouse proving ring and the combined 
nonlinearity and hysteresis was within ±0. 1 percent of the load cell full-scale reading. 


PROCEDURE 

A series of tests were made to determine the flow coefficient of the airflow mea- 
suring station, the flow and thrust coefficients of the four fixed -area nozzles, and ex- 
perimental terms in the gross thrust equation. 

The inlet (plenum) and altitude chamber pressures were varied to cover a wide 
range of pressure ratio (P /P^), while the inlet air temperature was set about 311 K 
(560° R) for all tests. (All symbols are defined in appendix A.) 


Airflow Calibration 

The flow rate calibration tests consist of two parts: one on the airflow measuring 
station, and the other on the four fixed-area nozzles. The inlet pressure P p was varied 
from 34. 5 to 144. 8 kilonewtons per square meter absolute (5 to 21 psia), and, at each 
P , the altitude chamber pressure was varied to give a range of P /P Q from 1 . 2 to 30 
as shown in table I . 


TABLE I. - AIRFLOW CALIBRATION 
TEST PRESSURES (NOMINAL) 


Plenum pressure 

Pressure ratio, 
P p /P 0 to P 7 /P Q 

kN/m 2 

psia 

34.5 

5.0 

1.2 to 15.0 

101.4 

14.7 

1.2 to 24.0 

124.1 

18.0 

1 . 2 to 27 

144.8 

21.0 

1 . 5 to 30 
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Airflow Measuring Station Calibration 


The accuracy of engine or component tests rests on the ability to determine accur- 
ately the airflow passing through the test article. One frequently used method is based 
on the continuity equation. For a compressible but perfect gas, the continuity equation 
can be expressed as 


W = AP 


2gy 

_RT(y - 1) 




( 1 ) 


At the throat of the bellmouth (station 1) , static and total pressures are recorded by 
probes located as shown in figure 4(c). There is one rake with five total pressure probes 
across the mainflow area of the throat, and there are four rakes each with six total pres- 
sure probes for boundary -layer jneasur ements plus four wall static pressure taps. 

An average total pressure P ^ at station 1 was obtained from the relation 

F i'^/ p i dA <2) 

The actual procedure in obtaining P^ was as follows: P^(l) was evaluated using equa- 
tion (2) with values of the five main flow area probes plus one set of values from one 
boundary layer rake . This calculation was repeated for each of the remaining three 
boundary -layer rakes. Then, P^ was equal to the average of these four values 

P^ = ijp^l) + P t (2) + P x (3) + P^)] (3) 


Figure 5 shows a typical plot of total pressure as a function of area. The last points on 
the curved portion are the boundary layer rake values. 

Using Pj as the total pressure in equation (1), the actual airflow W^, at the bell- 
mouth throat was calculated. Then using the measured total pressure P^ (not in- 
cluding the boundary layer values) at station 1 and equation (1), an ideal airflow, 
was calculated. An ideal airflow was also calculated using the average total pressure 
Pp in the plenum and equation (1). 
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0 20 40 60 80 100 120 140 

Area, in. 2 


Figure 5. - Station 1 total pressure including boundary-layer rake total pressures as function of 
area. 


In order to compare values of air flow under different inlet and exit conditions, a 
dimensionless flow parameter is used. It is known commonly as the discharge coeffi- 
cient C D . It is defined as 


C D,bm 


W 

Calculated actual air flow 1 
Calculated ideal air flow ^ 


( 4 ) 


Values of the discharge coefficient were calculated using ideal flow based on total 
pressure i( j at station 1 as shown in figure 6(a) . Values of the discharge coefficient 



Reynolds number index (station 1), REI 


(b) Based on plenum total pressure. 

Figure 6. - Airflow measuring station flow coefficient as function of 
Reynolds number index for cold-flow data 
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TABLE II. - FLOW RATE CALCULATION USING METHOD OF 


ANNULUS AREA SUMMATION 


Number of 
annulus , 
K 

Total pressure in annulus, 
P K 

Airflow in each annulus , 
W K 

kN/m 2 

psia 

kg /sec 

lb/sec 

1 

131 

828 

19. 

120 

1.0458 

2.3056 

2 
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15 

131 

.814 

19. 

118 

1.0454 

2.3047 

16 

131 

800 

19. 

116 

1.0450 

2.3038 

17 

131 

655 

19. 

095 

1.0407 

2. 2943 

18 

131 

276 

19. 

040 

1.0293 

2.2693 

19 

129 

966 

18. 

850 

.9891 

2. 1807 

20 

124 

519 

18. 

060 

.7999 

1.7634 


Annulus area, — (993. 16 cm) = (49.658 cm 2 ) (or 7.697 in. 2 ) 
20 

Static pressure station 1, p = 114.33 kN/m 2 (or 16. 5822 psi) 
Inlet air temperature, 310.69 K (or 559.24° R) 


Airflow in each annulus, 



Total airflow station 1, 


W ! = 


20 

W K = 20. 59 kg/sec (or 45.395 Ib/sec) 
k=l 
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TABLE III. - COMPARISON OF BELLMOUTH FLOW COEFFICIENTS 


Nozzle throat area, 
A 8 

Bellmouth throat 
static pressure, 
Pi 

Bellmouth flow coefficient, 
C D,bm 

cm 2 

in. 2 

kN/m 2 

psia 

Simplified 

method 

Summation 

method 

694 

107.5 

114.33 

16.58 

0.98662 

0.98433 

694 

107.5 

108.11 

15.68 

.98406 

.98785 

843 

130.7 

113.54 

16.47 

.98648 

.98878 

843 

130.7 

58.34 

8.89 

.98331 

.98269 

1073 

166.3 

100.92 

14.64 

.98420 

.98506 

1073 

166.3 

80.66 

11.70 

.98598 

.98451 

1186 

183.8 

74.95 

10.87 

.98445 

.98706 

1186 

183.8 

93.76 

13.60 

.98577 

.98888 


Average of eight headings 

0.98511 

0.98615 


were also calculated using values of plenum total pressure P p for ideal flow as shown 
in figure 6(b) . 

In an effort to check the accuracy of the actual airflow calculation using equation (1) , 
an actual airflow was calculated using the normal method of summation of annulus area. 
The flow area at station 1 was divided into 20 equal annulus areas. The air flow was 
calculated for each of these areas and summed to give the air flow for the whole duct. 
The results are shown in table II. A comparison of the results obtained using the two 
different methods is shown in table m. 


Exhaust Nozzle Calibration 

An accurate calibration of air flow through the exhaust nozzle is required in many 
engine programs. 

The procedure for determining the discharge coefficients for the four fixed-area 
nozzles is essentially the same as that for the inlet bellmouth. The actual air flow at 
the throat of the exhaust nozzle is equal to that calculated for the bellmouth if there is 
no leakage. Therefore, Wg = W^. 

Using values measured at station 7 for total pressure and temperature and insuring 
that the exhaust nozzle is always choked, the ideal air flow at the nozzle throat can be 
calculated using the following equation: 
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w„ Xl/r ' 1 


8, id 


RT ? \y + 1, 


, 1/2 


2gRT 


1 - 


1 y - 1 \ y + ly 


(5) 


The exhaust nozzle discharge coefficient is therefore defined as 


W, 


n w 
W 8, id 


(6) 


Thrust Equation Evaluation 


One of the important engine performance parameters is the gross thrust. In the 
engine test facility, the load cell in the test stand records overall thrust exerted on it by 
the entire system. The overall thrust hereafter called F , the measured thrust in- 
cludes, algebraically, (1) the gross thrust, (2) the pressure forces acting on the air in- 
let device, (3) the forces acting on the bulkhead seal due to pressure differential, and 
(4) the other miscellaneous facility system forces and restraints such as airflow leakage, 
structure joint friction, connecting piping, electrical lines, load cell hysteresis, and so 
forth. With proper facility design all the miscellaneous forces in item (4) can be made 
negligibly small, and are accounted for by a static calibration of the thrust measuring 
system. 

The measured thrust of the test system can be defined by the equation 


F = 
m 


M ( 


+ P 9 A 9 ' 


P 0 A 9 


(rav). + p. A. - PA. 
in m in p in 


A s< P p - P 0> + F 


fr.s tF d 


m 


Since the gross thrust is defined as 


and letting 


F = 
g 


(mv) 9 + p 9 A 9 - PqAq 


FIN = 


(mv) in + Pin A in “ P p A in 


(8) 


(9) 
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then 


F = F + FIN + A {P - P n ) + F„ + F, 
g m L s P O' fr,s d 


(10) 


Figure 7 shows the location and direction of these forces for the cold flow tests. A 
mechanical preload was used to insure a positive loading on the load cell at all times . 



Figure 7. - Thrust system showing location and direction of thrust terms used in gross thrust equation for cold-flow tests. 

TABLE IV. - THRUST CALIBRATION TEST 
PRESSURES (NOMINAL) 


Plenum pressure 

Pressure differential, 

P - P„ 

P o 

kN/m 2 

psia 

kN/m 2 

psi 

34.5 

5.0 

6.9 to 31 

1.0 to 4.5 

101.4 

14.7 

13.8 to 96.5 

2.0 to 14 

124.1 

18.0 

13.8 to 120 

2.0 to 17.5 

144.8 

21.0 

17.2 to 141 

2.5 to 20.5 
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The term FIN and the bracketed term in equation (10) must be determined in order 

to properly evaluate the gross thrust of the test article . 

The thrust equation evaluation program tests covered a range of inlet pressure P 

o P 

from 34. 5 to 144. 8 kN/m (5 to 21 psia), and of altitude chamber pressure Pq such that 

(P - Pq) varied from 6.9 to 138.0 kN/m 2 (1 to 20 psi) as shown in table IV. 

Bellmouth Pressure Forces 

The FIN term in the right side of equation (10) represents the momentum and the 
forces exerted on the bellmouth inlet due to the varying static pressure differential act- 
ing on the large and varying area curved surface. 

A theoretical determination of the combined bellmouth inlet pressure force term, 

FIN = (mv). + p. A. -PA. was made in terms of the corrected airflow rate at the 
' 'in in in p in 

bellmouth inlet. Therefore, 


— =A ^(Q-1) (11) 

P 

P 


where 


Q 


Gc 2 — 
g 




= yM: 


in 


in 


in 



( 12 ) 


The parameter Q calculated as a function of Gc and p in /P in is given in refer- 
ence 8. Thus, FIN can be calculated once the flow rates are known. 

Forces on the bellmouth seal and drag forces . - The second bracketed term in equa- 
tion (10) is composed of forces due to the difference in inlet plenum pressure and altitude 
chamber pressure and the restraining or frictional force produced by air flowing through 
the seal gap. There are also base drag forces included in this second term. The com- 
bined seal force plus drag force can be determined experimentally. To do this, the duct 
was blanked downstream of the airflow measuring station (station 1), and pressures in 
the plenum and altitude chamber were set as shown in table IV. Since there is no mass 
flowthrough the duct, (mv)g = 0 and P9 = Pg = Pq. Therefore, F^ = 0. 

Furthermore, (mv)^ n = 0 and p^ = P . Equation (4) is reduced to 

F m = -[VV P 0> + F fr,s + F d] < 13 > 


15 


That is the measured thrust in this test setup consists of the combined A g A p and fric- 
tion seal force plus drag forces. After FIN and the bracketed term are determined, 
the gross thrust can be determined. 

Gross thrust coefficient . - The ideal thrust equation is as follows: 

F id = C F,idW D> „ < 14 > 


where 


F, id 


r 

__ 

-> 

x y+l/y-1 

/n \ y ~ 1/V 


2 ( 2 \ 


> 

y - 1 v y + i / 

[ w J 






1/2 


p 9 ~ P 0 A 9 
' P 7 A 8 


(15) 


For any given nozzle and fixed inlet conditions of y, pg, P^, Ag, Ag, and Pg 


P n ^-q 

= (Constant) - — — - 

^ 9 PA 

7 8 


(16) 


Thus the ideal slope when plotted as a function of Pq/P^ is a straight line with a 
slope equal to the nozzle area ratio Ag/Ag. 

We can also define a similar coefficient with the expression 


F g ~ C F P 7 A 8 C D,n 


(17) 


or 




P 7 A 8 C D,n 


(18) 


where C p is the actual gross thrust coefficient. 

When the flow is expanded isentropically in the diverging section of the nozzle, the 
nozzle exit static pressure is fixed by the inlet total pressure and the nozzle area ratio 
If the altitude chamber pressure Pq is equal to or less than the nozzle exit pressure 
pg, the nozzle will flow full. As Pq increases above this value, flow separation will 
occur. When separation occurs the flow is apt to be more erratic and there is a loss 
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connected with the associated shock waves . The variation of C F against Pq/p,j should 
then reflect this and will no longer be a straight line. We have used this feature as an 
indication of nozzle flow separation, although direct observation was not made to con- 
firm it. 


Airflow Error Analysis 

Efforts were made to determine the resulting probable error in airflow computation 
due to various instrumentation errors. In general, the probable error is about ±0.4 
percent at the bellmouth throat design Mach number of 0.6. Details are given in appen- 
dix B. 


Gross Thrust Error Analysis 

Similarly an error analysis was made on the measurement of gross thrust. Calcu- 
lations indicated that the probable error was about ±0.3 percent at the load cell full- 
scale reading. Details are given in appendix C. 


Hot-Gas Flow Tests 

After the cold-flow calibration test runs were completed, a J85 turbojet engine was 
installed in the test stand. Hot-gas test runs were made with similar plenum and alti- 
tude chamber pressure settings up to exhaust gas temperatures of 1944 K (3500° R). 
Isentropic thrust coefficients and gross thrust coefficients were calculated for the ex- 
haust nozzles. 


RESULTS 

The results are presented in three parts: first, the results of the cold-flow cali- 
bration of the airflow measuring station, and the thrust and flow coefficients of the four 
fixed-area nozzles; second, the results of the hot-gas tests with a J85 engine which de- 
termined gross thrust and thrust coefficients for the four nozzles; and, third, a com- 
parison of the isentropic thrust coefficients for the hot -gas and cold-flow tests using the 
cold -flow data as a base line. 
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Cold-Flow Tests 


Flow measuring station (bellmouth inlet) calibration . - The average total pressure 
at station 1, P^ was evaluated by area integration of the recorded total pressure. A 
typical plot of against flow area is shown in figure 5. 

The actual flow at the airflow measuring station was calculated using P^ and the 
continuity equation. The actual airflow was also calculated by the method of summation 
of areas to check the method using the continuity equation and a single value of the total 
pressure. Table in compares several flow coefficients calculated using the airflow from 
the two methods. Agreement between them was generally within ±0. 2 percent. 

Airflow measuring station flow coefficients . - The bellmouth flow coefficients were 
computed by dividing the calculated actual airflow using P^ and p^ by the ideal flow 
rate, calculated using the free -stream total pressure without B/L probes at station 1, 
Pj id and pj. The results are presented in figure 6(a). 

The bellmouth flow coefficients were also calculated by dividing the actual airflow 
by the ideal flow rate calculated using the average plenum total pressure P p and p^ . 
The results are presented in figure 6(b) . 

Comparison of the flow coefficients obtained by the two methods showed that agree- 
ment was excellent. With either method the scattering of data points were all within a 
band of about ±0. 4 percent. The difference between the average values from the two 
methods is less than 0.07 percent. 

It was therefore concluded that using the plenum pressure for the calculation of the 
ideal air flow would be quite accurate. Thus, the total pressure rake at station 1 
could be eliminated for the hot-gas tests. This would provide a smoother airflow into 
the engine inlet as long as there was uniform bellmouth inlet conditions. 

Exhaust nozzle flow calibrations . - The four fixed -area water-cooled exhaust noz- 
zles are all convergent -divergent with an exit to throat area ratio (Ag/Ag) of about 1.3. 
The dimensions are given in figure 3(b) . These four nozzles were also used in the test 
program reported in reference 1. The pressure settings in the calibration runs are 
given in table I. 

The total air flow rate of station 8 was calculated with equation (5) using average 
total pressure at station 7) P 7 , id ‘ Figure 8 is a plot of typical recorded total pres- 
sures at nozzle inlet against the radii where the probes were located. These total pres- 
sures were integrated with respect to flow area and averaged to give P^ - d . In cali- 
bration runs, the actual air flow rate through the nozzle is assumed to be the same as 
the actual flow rate at station 1, that is Wg = W^. 
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Nozzle pressure 
ratio, 
P7/P0 


o L96 
□ 16. 15 



(a) Nozzle throat area, 843 square centimeters (130.7 in. 2). 



8642(£246 8 

Probe depth from centerline, in. 


(b) Nozzle throat area, 694 square centimeters (107.5 in . ^). 

Figure 8. - Typical nozzle inlet (station 7) total pressure profiles for cold-flow data. Plenum total pressure, 
atmospheric. 


The exhaust nozzle discharge coefficients were computed using the definition that 



These coefficients are presented in figure 9. The scattering of data points are within 
a ±0. 5 percent band except for the nozzle with a throat area of 694 square centimeters 
(107. 5 in. **) which is about ±0.7 percent. The average value of the discharge coefficient 
for each of the four nozzles varies from 0.9946 to 1.0025. Although it is quite common 
to have nozzle discharge coefficients greater than 1 for several valid reasons as dis- 
cussed in reference 9, while flowing hot gases, no clear cut reason can be given for 
these values greater than 1 determined using air at room temperature. However, data 
system accuracy and installation effects are the probable reasons for obtaining these 
values. 
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' 0 4 8 12 16 20 24 28 32 

Nozzle pressure ratio, P 7 /Pq 

(d) Nozzle throat area, 1186 square centimeters (183 8 in. 2 ). 

Figure 9. - Nozzle flow coefficient as function of nozzle pressure 
ratio for cold flow data. 


Thrust Measurements 

The exhaust nozzle thrust calibration and evaluation tests were conducted with the 
nominal pressure settings shown in table IV. The gross thrust is defined by equa- 
tion (10) . 

Combined seal force term plus base drag force . - Physically the pressure probes 
for P were located in the plenum slightly ahead of the labyrinth seal, and probes for 

r 

Pq were in the altitude chamber at a location roughly where the throat of the exhaust 
nozzle would have been as shown in figure 2(a). 

A typical result is shown in figure 10 where the experimentally determined seal 
force plus drag force is plotted against (P - Pq). It is seen that the data follow an al- 
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Seal pressure differential (P p - Pq), psi 

Figure 10. - Combined seal force and drag force as function of pressure 
differential. 



most straight line which passes through the origin. This indicates generally good seal 
design and linearity of the system. Also shown on this figure is the calculated A g AP 
based on the physical seal area (A g ) and the measured AP(P p - P^). Figure 11 shows 
the difference between the measured and calculated seal force (friction) and base drag 
force term plotted against (P - Pq) . This term is obtained by subtracting the calcu- 
lated seal force from the experimentally measured force and is an indication of the drag 
component included in the overall measured term. 


Bellmouth Pressure Forces 

The combined forces acting on the curved and varying-area surface of the bell- 
mouth was calculated by equation (11). Figure 12 shows the combined bellmouth inlet 

o 

forces divided by the plenum pressure as a function of the corrected airflow. 
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Figure 12. - Theoretical and experimental combined inlet force over plenum total pressure as function of 
corrected inlet airflow. 


Gross Thrust 

The sum of the measured thrust, the seal force plus drag force, and the combined 
bellmouth inlet force determines the gross thrust produced. A typical plot of gross 
thrust against exhaust nozzle pressure ratio, P^/Pq, with inlet pressure constant for 
cold flow is presented in figure 13. 
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Figure 13. - Gross thrust (from cold flow data) as function of nozzle 
pressure ratio. Nozzle throat area. 694 square centimeters (107.5 
in. 2 ). Inlet pressure held constant 


Isentropic Thrust Coefficient 

The ratio of gross thrust to the maximum thrust produced under isentropic condi- 
tions is defined as the isentropic thrust coefficient and is a measure of the nozzle thrust 
efficiency. These parameters for the range of nozzle pressure ratios tested are pre- 
sented in figure 14. 


Actual Gross Thrust Coefficient 

The commonly defined thrust coefficient, Cp ^ in the ideal thrust equation 

F id = C F,id P 7 A 8 C D,n ( 20 ) 

is a straight line -with a slope of Ag/Ag when plotted as a function of Pq/P^. The cal- 
culated gross thrust can define a similar parameter called the actual gross thrust coef- 
ficient. A typical comparison of these two coefficients is shown in figure 15. At low 
values of Pq/P^, the C F values cluster quite close to the slope line of the ideal case 
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and this serves as a check that the nozzle is flowing full . Deviation of the gross thrust 
coefficient values from the ideal slope line where Pq/P^ ~ 0.35 is an indication that 
the exhaust gas is separated from the wall and the nozzle is no longer flowing full. This 
detached flow can be highly asymmetrical as well as unsteady. Data beyond this point are 
probably not applicable in a steady-state performance evaluation program. 


Hot-Gas Flow Tests 

In the hot-gas test runs, a J85 turbojet engine was installed as a gas generator, and 

tests were run with exhaust gas temperatures up to 1944 K (3500° F). For the 694 and 

843 square centimeters (107. 5 and 130. 7 in. nozzles, the temperatures at the nozzle 

inlet were measured with temperature rakes. For the 1073 and 1186 square centimeters 

2 

(166.3 and 183. 8 in. ) nozzles, the temperatures were calculated using choked flow re- 
lations due to difficulty of measuring high temperatures resulting from afterburner op- 
eration. 

The gross thrust was again the sum of the measured thrust, the combined bell- 
mouth inlet force, and the seal force plus drag force. Using gross thrust values, the 
isentropic thrust coefficients were calculated and are presented in figure 16 for each of 
the nozzles. 

Actual gross thrust coefficients were calculated and a typical result is presented in 
figure 17 for the 694-square-centimeter (107. 5-in. nozzle. As in the case of the cold- 
flow data, the hot-gas data follows the slope of the ideal case at low values of Pq/P^ . 
Deviation of the gross thrust coefficient values from the ideal slope at values of 
Pq/P^ rj 0.35 is again an indication of nozzle flow separation. 


Comparison of Hot-Gas and Cold-Flow Results 

One of the prime purposes of these tests was to determine how accurately gross 

thrust could be evaluated for this turbojet engine and to determine if the method used 

for measuring thrust in these tests could also be used for other engine and nozzle tests. 

Experimental nozzle tests reported in references 9 and 10 showed that at a higher 

nozzle inlet temperature, the nozzle flow coefficient value was less than the value 

determined from the cold flow data. Because of these results an attempt was made to 

calculate hot-gas flow coefficients for the fixed-area nozzles from the hot-gas data. 

Figure 18 presents a typical curve of hot-gas flow coefficients as a function of nozzle 

2 

pressure ratio for the 694-square-centimeter (107. 5-in. ) nozzle. These flow coeffi- 
cients were used in determining the thrust coefficients for the four nozzles for the hot- 


25 



After burning 
hot gas data 
Nonafterburn mg 
hot gas data 



5 

’ O tA 


Q. C 
. CD 
CD O 


c 5 
OJ S 


S*s 

s e 

CT3 .C 
> Q_ ■«— 

S? 

-a. M 

5 M 



26 


Nozzle pressure ratio, P 7 /Pq 

(dl Nozzle throat area, 1186 square centimeters (183.8 in. 2 ). 

Figure 16. - Nozzle isentropic thrust coefficients for hot-gas 
flow with exhaust gas temperatures up to 1944 K, 350CP R 
as function of nozzle pressure ratio. 


Pressure ratio at Temperature ratio at 
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P 7 /P 7 T 7 /T 7 



Figure 18. - Nozzle flow coefficients for 694-square centimeters (107.5- 
in.2) nozzle from hot-gas data as function of nozzle pressure ratio. 
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Figure 19. - Temperature ratios and total pressure ratios at station 7 as function of location 
at station 7. Temperatures measured with water-cooled traversing probe for afterburner 
operation. 
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gas tests. Also, figure 19 shows that the temperature profile at the nozzle inlet for hot- 
gas flow up to 1944 K (3500° R is quite different from the temperature profile for the 
cold-flow tests with ambient temperature at the nozzle inlet. The total pressure pro- 
files at the nozzle inlet for hot-gas and cold-flow operation are presented in figure 19. 

The isentropic thrust coefficients for the hot-gas and cold-flow tests as a function 
of nozzle pressure ratio are presented in figure 20. There is favorable agreement be- 
tween these values for the hot-gas and for cold-flow data. For all the nozzles on the 
average, the hot-gas coefficients were within ±0. 59 percent of the cold-flow base-line 
data. 
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Hot-gas data 
Cold (low base line 
Plain symbols denote 
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(a) Nozzle throat area, 694 square centimeters (107.5 in . \ 





(d) Nozzle throat area, 1186 square centimeters (183. 8 in.^l. 


Figure 20. - Comparison of nozzle isentropic thrust coeffi- 
cients for hot-gas flow and cold flow as function of nozzle 
pressure ratio 
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SUMMARY OF RESULTS 


The cold-flow tests and the hot-gas tests conducted in this program have yielded 
the following principle conclusions: 

1. In a well-designed test setup with an accurately calibrated thrust measuring sys- 
tem, gross thrust can be determined within 0. 6 percent of the load cell full-scale read- 
ing for both hot -gas and cold-flow operation based on a probable error analysis. 

2. A direct comparison of the nozzle isentropic thrust coefficients for both the hot- 
gas and cold-flow tests using the cold-flow data as a base line showed that agreement 
between these coefficients was on the average within ±0. 59 percent for the four nozzles. 
The results indicate that the method used to determine gross thrust for the turbojet en- 
gine in this test could also be used in a future engine and nozzle evaluation tests. 

3. In the calabration of the airflow measuring station, the bellmouth flow coefficient 
was calculated by two different methods. In one, the free -stream total pressure outside 
the boundary layer at station 1 was used; in the other, the average total plenum pres- 
sure was used. The results from the two methods are within 0.07 percent of each other. 
It is therefore permissible to eliminate total pressure probes at the bellmouth throat 
(station 1) during hot -gas runs for engine testing, when uniform flow conditions are pre- 
sent at the bellmouth inlet. This can help improve airflow pattern into the engine by 
eliminating possible wakes from the instrument rakes. 

4. The flow coefficient for the airflow measuring station (station 1) was 0.984 with 
a scattering of ±0.4 percent. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 30, 1971, 

720-03. 


29 



APPENDIX A 


SYMBOLS 


A area 

a length (nozzle inlet to throat) 

abs absolute 

b length (nozzle exit to throat) 

C coefficient 

D diameter 

F thrust or force 

Gc mass flow function 

g gravitation constant 

M Mach number 

mv momentum 

P total pressure 

p static pressure 

Q stream thrust function 

R gas constant 

REI Reynolds number index 6 /<p V(T 

T total temperature 

t static temperature 

V velocity 

W flow rate 

z compressibility factor 

y ratio of specific heats 

2 

6 ratio of total pressure to NASA standard sea-level pressure of 101 kN/m 
(14.696 psia) 

6 ratio of total temperature to NASA standard sea-level temperature of 288 K 
(518.67° R) 

p density 
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(p ratio of viscosity to the viscosity at NASA standard sea-level conditions 

Subscripts: 

a A* afterburner 

bm bellmouth 

D flow or discharge 

d drag 

exit nozzle exit, station 9 

F thrust 

fr friction 

g gross 

id ideal (one dimensional) 

in inlet station of bellmouth 

K annulus number 

L local 

m measured 

n exhaust nozzle 

P plenum 

s seal 

si sea-level 

0 altitude chamber 

1 bellmouth throat 

7 exhaust nozzle inlet 

8 exhaust nozzle throat 

9 exhaust nozzle exit 
Superscripts: 

average 

* critical state (M = 1) 


APPENDIX B 


DETERMINATION OF AIRFLOW ERROR 

A calculation was made to determine the error in airflow due to assuming y = 1.4. 
For the design Mach of 0. 6 for the bellmouth throat and for the range of temperatures 
covered, the deviation was approximately +0.028 percent. 

The cumulative error due to various measurement and recording errors from indi- 
vidual instrumentation is in table V. The error given for each individual item was sup- 
plied by the manfacturer . 


TABLE V. - INDIVIDUAL AND CUMULATIVE 
ERRORS OF INSTRUMENTATION 


[Values in percent. ] 



Error 

Cumulative error 

Damper 

0 . 1 

0 . 1 

Load cell 

0. 1 

1 0. 105 

Digital counter 

0.1 

/ 

Thermocouple 

0. 375 


Oven 

0.2 

)• 0. 45 

Average digital error 

0. 15 

J 


The probable error resulting in the airflow calculation due to instrumentation errors 
shown in table V can be estimated as follows: 


W 1 ' C d,bm W l, id “ C d,bm A l P l,id 


2 gy 

RT( r - 1)1 


! 1/2 



x 



(Bl) 
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Letting 



I 

Differentiating equation (B2), then dividing by and letting (dy/y) = (Ay/y) , yield 
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-1 


r 



For bellmouth airflow Mach number = °- 6 ’ P l,id = 101.4 kilonewtons per square 
meter absolute (14.7 psia) p^ = 79.3 kilonewtons per square meter absolute (11. 5 psia) 


Pi 


P 


1, id 


0.78 


AC , , 

— i_bm _ 0.005 

^d,bm 


Ap t 

Pi 


0.001 


1, id 
l,id 


0.001 
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= 0.001 (assumed) 


i= 0.0045 

T 1 

AP 1 ^ = 0 -° 01p i id = 0.1014 kN/m 

A Pl = O.OOlPj - 0.0793 kN/m 2 (i 

Substituting these values and constants k 0 and k„ 
o 2 to 

(ep) = (0. 00610) or a probable error of 0. 61 percent. 


(0.0147 psi) 

.0115 psi) 

into equation (B4) gives 
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APPENDIX C 


RESULTANT ERROR IN GROSS THRUST MEASUREMENTS DUE TO 
INSTRUMENTATION ERRORS 


Using equation (10), the gross thrust is 

F = F + Rmv). + P. A. - P A. 1 + [a (P - P n ) + F. + F ,~1 
g m 'in in in p inj L s' p 0' fr,s dj 

Since the quantity [ A s< P p - P 0> + F fr,s + = F g is measured as one quantity 

(see fig. 10), it is not necessary to evaluate the errors contributed by A P P n , 

s p u 

Ffr g , and F^ separately. Instead, the measurement error in F g as a whole will be 
used. 

Therefore, letting £(mv) in + P in A in - P p A. J = FIN, 


F = F + F + FIN 
g m s 

Substituting the definition of FIN (eq. (11)) yields 

F g = F m + F s + P p A m<« - » 

The absolute value of the error in gross thrust measurement is 

AF = AF + AF + P AA. + A. AP 
g m s p in in p 

The probable error expressed as a fraction of the gross thrust measurement is 
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Assume P = 144.8xlO J newtons per square meter absolute (21 psia), = 22 240N 
(5000 1b), A in = 0.6819 square meter (1057 in. 2 ), M = 0. 6, AF = ±22. 24N(±5 lb), 

AP = 144.8 newtons per square meter absolute (0.021 psia), and AA = 0.6819X10 2 


square meter (1.057 in. 2 ). Substituting these values into equation (Cl) results 


in 



— [(5.0) 2 

(5000) 2 L 


+ (5.0) 2 + (1057x0. 021) 2 


+ (1.057x21) 


2 ' 



I X1036 

(5000) 2 



0.000041 


(ep) 2 


(0. 006 4 ) 2 


or a probable error of 0. 64 percent. 
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